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Abstract 



We calculate limits to the properties of massive, unstable neutrinos using 
data from gamma-ray detectors on the Pioneer Venus Orbiter Satellite; a 
massive neutrino emitted from SN1987A that decayed in flight and produced 
gamma rays would be detectable by this instruments. The lack of such a signal 
allows us to constrain the branching ratio to photons (By), mass (m u ), and 
radiative lifetime (r 7 = t/B^). For low mass (m < T ~ 8 MeV) neutrinos 
decaying v — ► 1/7, B 7 < 3 x 10 -7 , for m u r < 10 6 keV sec, and 6 7 < 6 x 
10~ 14 m u T/ keV sec for m v T > 10 6 keV sec; limits for high-mass neutrinos 
are somewhat weaker due to Boltzmann suppression. We also calculate limits 
for decays that produce gamma rays through the bremsstrahlung channel, 
v — > z/e + e~7- In the case that neutrino mass states are nearly degenerate, 
5m 2 /m 2 <C 1, our limits for the mode v — ► 1/7 become more stringent by 
a factor of 8m 2 /to 2 , because more of the decay photons are shifted into the 
PVO detector energy window. 
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I. INTRODUCTION 



The occurrence of Supernova 1987A in the Large Magellanic Cloud has proven to be 
among the most fruitful experiments in the heavenly laboratory for confirming "known" 
physics and constraining new physics. Aside from its obvious impact upon the study of 
the late stages of stellar evolution in general and upon supernova physics in particular, 
models for SN1987A have become a test bed for the study of the couplings of light particles 
(neutrinos, axions) to ordinary matter |TJ. In this work, we discuss limits on the properties 
of neutrinos independent of a specific model for the supernova based upon the emission of 
thermal neutrinos from the hot nascent neutron star. 

When a supernova occurs, the bulk of the binding energy of the progenitor star (~ 
3 x 10 53 erg) is released in neutrinos, as predicted by theory and confirmed by the observation 
of a neutrino burst from SN1987A, with a characteristic temperature of about T v « 4.5 MeV 
(for low-mass electron neutrinos; other species are predicted to have a higher temperature, 
T ~ 8 MeV, because of their different coupling to the prevalent electrons J^fl)- If & t least 
one species of neutrinos is unstable and couples to the photon, then some of these neutrinos 
will decay to photons en route, which are potentially detectable as MeV gamma rays. At 
the time of the supernova burst's arrival at earth and environs, there were several satellites 
operating in the solar system capable of detecting the decay photons in the course of their 
watch for gamma-ray bursts. Analyses of the data from one of these detectors, on board the 
Solar Max Mission (SMM) satellite, has already been presented @. Here, we examine the 
data from the Gamma Burst Detector on the Pioneer Venus Orbiter (PVO). While the PVO 
detector was smaller, and its energy window is not well-matched to that of the supernova 
neutrinos, it had 47r acceptance and was in an environment free of the Earth's radiation belts 
(leading to lower backgrounds). In addition, more high quality data is available (> 8000 sec 
vs. 10 sec for SMM). 

The paper is organized as follows. In the next section an exact formula for the expected 
gamma-ray flux is derived and important approximations are developed. In Section |TJ, the 
PVO data are discussed and rigorous limits are derived in the simplest regime. The next 
four sections build upon these results, expanding to more complicated regimes. The final 
section is a brief summary. 



II. GAMMA-RAY SIGNAL 



We can write the expected fluence of gamma rays from decaying neutrinos with mass m u 
and mean lifetime r as || 



dN 



7 * v ' dE n(fj)dfj. - 



-td/y 



x6\t — t d 



-dt n 



l- Vf i+^-\y/l-(vt d /DY(l 



dt. 



(1) 



where is the fraction of decays that produce a gamma ray. The first factor is the overall 
flux of neutrinos from a supernova at a distance D. L#(E) is the differential number flux of 
neutrinos of energy E, so E T = J dE EL# is the total luminosity in neutrinos. The second 
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factor gives the fraction that decay into a "lab-frame" angle arccos \i. The third factor gives 
the fraction that decay at time td- Finally, the delta function selects the photons with a 
given td, E, /i that arrive at a time t after the arrival of massless neutrinos at the detector. 
The Lorentz factor is 7 = E/m u and the speed v = y/1 — 7~ 2 . 

The function n(p) depends on the distribution of daughter photons in the neutrino rest 
frame, and therefore on the particular decay channel involved. First, we consider the two- 
body decay, v — ► u'j. Because the neutrino is a spin-1/2 particle, and the photon a spin-1 
particle, this reaction can proceed in one of two ways: with the helicity of the daughter 
neutrino either parallel or antiparallel the initial helicity. From quantum mechanics, then, 
the distribution of the photon in the rest frame of the parent will be proportional to either 
(1 ± pj)/2, where p, is the cosine of the rest-frame angle between the directions of the parent 
neutrino and the photon. Transforming into the lab frame gives us the distribution 
Note that we have assumed the neutrinos are ejected from the supernova in an instantaneous 
burst; as long as the actual duration of the pulse is small compared to the timing resolution 
of the detector, which is the case, this is a good approximation. 

Because we are not interested in the decay angle, but rather the photon energy, we write 
n(/j,)d[A = f(E, k)dk, where k is the gamma-ray energy, related to the decay angle by 

This gives 

f{E,k) = -^{EvTE±<2k), (3) 

for each of the helicity possibilities. (For reference, an isotropic decay would give f(E,k) = 
l/(Ev) = 1/p, where p is the neutrino momentum.) For low-mass neutrinos with d«1, 

( 2k/ E 2 no flip 

f(E,k) = l2(E-k)/E 2 flip . (4) 
[ 1/E isotropic 

Each of these should be multiplied by a Heaviside function Q(E — k) to require that the 
daughter photon be less energetic than the parent. Further, we require that the decay does 
not occur inside of the progenitor envelope which would considerably alter the energetics of 
the explosion and lead to an independent constraint which is important for short lifetimes 

The factor inside the delta function is especially complicated. This occurs because at 
any given time, the detector is receiving photons from neutrinos that have decayed in a 
complicated shape, approximately an ellipsoid with the supernova at one focus and the 
detector at another, further complicated by the speed v < 1 of the massive neutrinos. This 
includes photons that have left the supernova pointing far away from the detector but have 
decayed at large angles toward the detector. Obviously, for low mass neutrinos which leave 
the supernova at highly relativistic speeds, the fraction which take such a path is very small. 
To simplify this expression, we shall require that td "C D ~ 5 x 10 12 sec — most of the 
neutrinos decay well before they reach the earth. Otherwise, the flux is greatly reduced and 
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the limits are correspondingly weaker; the neutrinos may typically take a path that causes 
them to decay well after our observations end. In Sec. [VI I| we discuss long lifetimes. 

In this limit, the delta function becomes simply 8 (t — t<j(l — u/i)) = (2 r jk/m)5(t d — 
2jkt/m), and we can perform the integral over t d : 

dN = B "i L *(f) HE, k)—e- 2kt l m » T dt dE dk (5) 
4ttD 2 v 'm u T v ' 

(Similar expressions have also been derived in J|||.) We shall assume that the neutrino- 
number luminosity is given by a zero-chemical-potential Fermi-Dirac distribution with known 
temperature and total neutrino energy, a reasonable approximation ||. For now, we will 
present results for low- mass neutrinos {i.e., m u <C T v ), 

. 120 E T E 2 
L^E, mu = 0) = — — TT - m -, (6) 

where Ex — 10 53 erg is the total energy in one species of massless neutrinos. We treat the 
case m u > T u in Sec. [V] below. 

Finally, we can integrate the above expression over all neutrino energies E and over one 
time bin from t to t + St to get an expression for the spectrum of photons incident on the 
detector during that time interval. For m v <C T u , 

In this expression, the function h(k/T v ) results from the integral over the neutrino energies. 
It is of order unity for the paramater ranges of interest, and it is largest in the case of "no 
flip," which we will assume from now on since it gives the most conservative estimates of 
the parameters. In that case, it is given by h(y) = yln(l + e~ y ). 

Although this signal depends nonlinearly on the parameter m u T, the expression simplifies 
when m u r is much greater or less than kt or kSt, where A; is a typical photon energy, giving 

6(k t) - -^L- — ^h(klT ) x { m ' m » T « k5t (8) 

mt) ~ AnD 2 7^Tr {k/K) \2k5t/m u r, m u r » kt ' W 

In the former case of small m u T, there is no appreciable relativistic delay before the decay 
of the neutrinos, so essentially all of the daughter photons arrive in the first time bin after 
the supernova. In the case of large m u T, the flux is essentially constant over the time of the 
observations, so the signal is proportional to the width of the time bin. Note that only in 
the latter case does the fluence actually depend on the value of m u r. 

In order to calculate the expected signal from the theoretical spectrum, we must fold 
that spectrum with the appropriate response function. The signal expecied in the zth energy 
channel is 

Si(t) = J dkRi(k)<l>(k,t) = J2^Mt) (9) 

3 

where Rij is the response of detector % in energy bin kj, and </>j(t) is the theoretical spec- 
trum averaged over energy bin j at time (or time bin) t. This, combined with a separate 
observation of the background rate in each of the four detectors, provides the theoretical 
signal to be compared with the observations. 
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III. GAMMA-RAY DATA AND ANALYSIS 



To obtain our limits we use data from the supernova with data from Pioneer Venus 
Orbiter Gamma Burst Detector (PVO GBD) ||, luckily in operation in February, 1987. 
The GBD has four energy channels, roughly 100 - 200 keV, 200 - 400 keV, 500 - 1000 keV, 
and 1 — 3 MeV towards the direction of the supernova, which was propitiously directly 
overhead at the time, giving the maximum effective area. We have data for about 1500 sec 
prior to the supernova lights arrival at Venus (for calculating the background), and for 8000 
sec after, for time bins of either 12 or 16 sec in duration. We show the data in Figure |l|. 
We have verified that there is no clear signal in any of the four channels; further, the data 
is consistent with a constant Poisson rate in each detector. 

To calculate limits on our parameters B 1 and m u r, we use the folded spectrum Si(t) 
and our measurement of the background rate (observed for a length of time tb) in each 
detector to construct a likelihood function given the observed gamma-ray counts in each 
detector. We assume that the folded spectrum Si(t) gives the mean of a Poisson process 
governing the detected number of counts, and that the rates for each detector are high 
enough to approximate this by an appropriate Normal distribution, for ease of calculation 
(in the 1 — 3 MeV bin, with the lowest fluence, there are approximately 40 counts per bin). 
This gives 



£(6) = J] N(D ir , kStj + S tj (6), a%) 



(10) 



where 



N(x; /i, <r 2 



V2 



: exp 



1 (x — n) A 



2 a 2 



(11) 



gives the Normal distribution, 6j is the background rate in detector i (observed for a time 
U), Stj is the length of time bin j, and Dy, Sij are, respectively, the observed and theoretical 
signal in those time bins, where the latter is calculated with the set of parameters represented 
by 9. Finally, the variance is given by afj = Sij + bidtj(l + Stj/tb), the sum of the theoretical 
variance of the signal and that due to the background rate (including a small contribution 
reflecting the uncertainty in that rate; this latter effect is somewhat more difficult to include 
if a Poisson distribution is explicitly used but is in any case negligeable) . We define a x 2 
statistic, 



X 



-2 In C + const = 



l] L 



In a; 



$ij n (bi8tj — Dij) Sij (Di 
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\2l 



(12) 



Note that the model is nonlinear, and the variance cr?- depends on the model parameters, 
so we have defined this quantity including the In a^- term and that the usual \ 2 distribution 
does not apply; instead we must apply Bayes' theorem and integrate over the likelihood to 
determine confidence intervals on our parameters. 

Because of the two terms contributing to the variance, the form of \ 2 depends on which 
dominates. For S%j ^> hStj, afj ~ Sij, and the Sfj/afj term dominates, so x 2 ~ Sij. When 
the neutrino signal is small, the background contribution dominates, and x 2 — const. These 
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regimes are shown in Figure |2], where we plot x 2 as a function of m v T for several values of 
B T We have used a neutrino temperature of 8 MeV, appropriate for mu and tau neutrinos. 

Immediately, we see the character of the limits on the parameters. For m u r < 10 7 keV sec, 
X 2 oc -B 7 ; in this regime only the data from the first time bin after the supernova contributes. 
Then x 2 oc B^/m u T for m v r / B^ < 10 13 keV sec; now, the full data set provides information. 
Finally, x 2 — const for m^r/B^ 10 13 keV sec; in this regime the background dominates over 
the theoretical signal. Note that this latter area of parameter space provides the "maximum 
likelihood" (or x 2 minimum); there is no neutrino signal and we calculate only limits on 
parameters. In fact, there is a slight deficit of counts with respect to the background 
calculated from the time before the supernova; otherwise we might expect to see a weak 
maximum likelihood somewhere in the large-m^r regime. 

In Figure [J[ we show a single contour of x 2 m the m v T-B^ plain. From a "Bayesian" 
viewpoint, this figure represents the the complete inference from the data; we obviously 
would prefer to quote limits on the parameter space. 

Because x 2 is constant for both large and small values of the parameters m u T, we must be 
careful in normalizing our probabilities. Without any other information to guide us, we might 
expect the appropriate prior to use would be either p{6) oc const or p{9) oc 1/9 corresponding 
to constant probability per linear and logarithmic interval, respectively; however, neither of 
these choices converge when integrated out to infinity multiplied by a constant likelihood. 
These priors, when normalized, assign a vanishing weight to any finite region of parameter 
space; usually, x 2 g° es to infinity like the square of some parameter in a linear gaussian 
model with independent errors, so the wings are automatically disfavored regardless of the 
prior. 

To get meaningful results in this case we must make sure that the interesting regime 
of —6 < log-B 7 < 0, 5 < \og(m u r/ keV sec) < 13 is given finite prior weight. This 
corresponds to choosing a prior with some cutoff outside of this region, and will result in 
a limit corresponding to a value of x 2 somewhere along the slope between its small and 
large values. To connect with other analyses, we will adopt the following "frequentist" 
procedure. We will assume that x 2 is distributed in a x 2 distribution considering the data 
as a frequentist random variable. We have 550 time bins and four detectors, so there 
are 2200 degrees of freedom. For this distribution, a one-sigma fluctuation corresponds to 
Ax 2 = 2230, a three-sigma (or 99%) fluctuation to Ax 2 = 2357; we choose the latter as our 
limit; from the shape of the likelihood function it is clear that any comparable Ax 2 will give 
similar bounds. We also note that a signal in the small-m^r regime may not be detectable 
with this algorithm; the absence of a local minimum in that region, however, implies that 
this should not be a significant worry. (We have also used a x 2 defined rigorously from the 
true Poisson distribution, which effects the limits on the parameters by somewhat under one 
order of magnitude throughout.) The allowed region is shown in Figure ^. It corresponds to 

B 7 <3x 10~ 7 m u r < 10 6 keV sec (13) 

77? 7~ 

K<2x 10~ 13 - m u r > 10 6 keV sec. (14) 

keV sec 

for neutrinos with a temperature of 8 MeV. (The limits scale roughly as T~ 2 .) The latter 
bound corresponds to m u T^ > 5 x 10 12 keV sec. This is less restrictive than the limits of 
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Oberauer et al. j|, due to the fact that the PVO GBD could only detect gamma rays with 
energies below 3 MeV, compared to 25 MeV for the SMM satellite. However, we believe this 
analysis to be more rigorous and the PVO data to be of higher quality. 



IV. OTHER DECAY MODES 

So far we have considered only photons produced from the simplest radiative decay of 
a massive neutrino species: v — > 1/7, but many other channels are possible. For these 
reactions, where the rest-frame photon energy is no longer given by the simple k = m v /2, 
we must allow for a distribution of decay products: f(k, p)dkdp, gives the fraction of photons 
produced with rest-frame energy k into angle /2 = cos 6?. Then, the final spectrum is given 
by 



dN 



1 B~ 



dkdt AttD 2 t 



dEL # (E) dfif 



k 



7(1 + Vfl) 



3T) A* 



-7(1+«m)*/ t 



(15) 



where we have still assumed that the decays occur near the supernova, and the integral is 
taken from k to 00. (We discuss the case of large masses, m v > T u in Sec. [V] below.) For 
the appropriate choice of f(k,fl) oc S(k — m u /2)fi(jl), we reproduce the earlier 2-body decay 
formula, Eq. (|7|). 

In particular, we consider the bremsstrahlung process, v — > z/e + e~7, where typically 
v — i/ T , v' — v e . Because this is no longer a two-body decay, the spectrum in angle and 
energy of the daughter photons is considerably more complicated, and there has been no 
exact calculation performed. Following Oberauer et al. |4j], we make several simplifying 
assumptions. First, we assume isotropy of the photons in the rest frame of the neutrinos; 
this is reasonable if the helicity states of the parent neutrinos are produced in equal num- 
bers. This gives f(k,p) = f{k)/2 (which still implicitly depends on p, through the Lorentz 
transformation to the lab frame). 

Up to factors of order unity, we assume after Oberauer et al. that the bremsstrahlung 
energy spectrum is given by 



dT 



br 



dk 



7T k 



a 1 
7r kr P 



(16) 



where r and r e refer to the process without a daughter photon: v u'e + e , allowing us to 
absorb a branching ratio factor into r e = r/B e . Now, we can write the gamma-ray flux as 



dN 
dkdt 



1 



a 



a 



AnD 2 m u T e 7r 



-(l+v)Et/m„T 
u=—(l—v)Et/m v T 



(17) 



If we assume that jt <C t, the angular average becomes unity. (A Maxwell-Boltzmann 
distribution for L#(E) then reproduces the expression of Eqs. (6-7) of Oberauer et al. [|J) 
With the additional assumption that a negligible fraction of the decays occur inside of the 



7 



progenitor and that detected gamma rays have energies k C T (i.e., E m[n < T), we get the 
simple formula 

dI^ = _E^J_al 
dkdt AtiD 2 m u T e it k 

where Et again gives the total energy in the decaying neutrino species. Note that this 
expression is correct for any form of L#(E). At early times the flux is constant, although 
there will be an exponential decay for yt ^ r. In the two-body case, the timescale for 
gamma-ray detection is set by m u r/k; in the bremsstrahlung case, by m u r/E ~ m u r/T. 

Assuming that most neutrinos will have E ~ T, the assumptions we have made require 
that t/r <C m v /T <C 1 for this formula to hold. Note that this is a restriction on the param- 
eter space of m v and r, since we have approximately 8000 sec of data, and the temperature 
is of order 8 MeV for any neutrino species. 

Because the flux is proportional to l/m u T e , we will be able to put limits on the combi- 
nation m v T e = m u T/B e . Moreover, because the dependence is the same as the large-m^r 
regime of the 2-body decay case (cf. Eq. (|])), the probability density will be of the same 
form. When the variance is dominated by the signal, x 2 ~ S{j ^ ^/ m u T e', when it is 
dominated by the background, x 2 ~ const (the same constant as in the 2-body case). The 
crossover occurs at m u T e « 10 15 keV sec. The value of x 2 f° r this model is shown in Figure |j. 
Using the same quasi-frequentist definition of a 99% confidence level gives limits of 

TD 

m v r e > 1.5 x 10 12 keV sec or < 7 x 10~ 13 keV -1 sec" 1 . (19) 

Because the bremsstrahlung spectrum peaks at a lower energy, and due to the long time 
baseline of the PVO data, this limit is comparable to other SN1987A limits for this decay 
channel and we believe more reliable, due to the higher-quality data set. 



V. VERY MASSIVE NEUTRINOS 

All of these expressions are considerable more complicated in the case m ^ T. We 
will still assume a zero-chemical-potential FD distribution, this time applying to massive 
particles: 



120 E^ EJE 2 + ml 

L # (E,m u ) = —fj [mjTj) ^ + qE/t / , (20) 

with a "suppression factor," j(x), along with the requirement that E > m u . The factor 
j(x) is just the usual Boltzmann suppression (j(x) oc x 3l/2 e~ x , for x 1); we use an 
approximation that is good for m/T < few, j(x) ~ exp (— 0.15a; 2 ). (Sigl & Turner [|?| have 
calculated the effect of the changing neutrinosphere temperature and radius on this naive 
expectation; the effect is small for m u < 40 MeV, at least for r > 10 -2 sec.) For low-mass 
neutrinos we assumed k > m u ; now, we can only integrate over neutrino energies greater 
than max(m^, k). In this expression, Et — 10 53 erg remains the total energy for the low-mass 



8 



case; the total energy released is E T j{m u /T u )\ which for large masses is less than 10 53 erg 
since j < 1 

In addition to the mass threshold effects, with massive neutrinos (m„ ^ k) we must now 
take into account the loss of any photons produced inside the envelope of the supernova, 
i?env = 100c sec. Thus, we require that vta > R env , or 



E>E P 



m v \ I + 



R P 



2kt 



(21) 



Note that E cnv > m u , so this supersedes the requirement that E > m u , but the re- 
quirement that E > k remains. Thus, we must integrate over neutrino energies from 
E min = max(fc, E env ). This integral, the equivalent of h(k/T) above, cannot be done in 
closed form, but again it can be approximated by a gaussian (at least for the isotropic case 
f(E,k) = l/p): 



dN 



B- 



dk dt AnD 2 



j (m u /T v ) 



120 Et 2k 
7n 4 T?, m,,T 



-2kt/m I/ T 



g(E miQ /T) 



with 



9ix) = f 



dy 



x 



l + e 1 



IT 

12' 



-0.2a; 2 



(22) 



(23) 



the factor of 0.2 in the exponent approximates the shape of the integral for x ^ few. This 
differs from the massless case by a total suppression factor 



j{m u /T u )g(E min /T) 
h(k/T) 



cxp 



-0.15 



m. 



0.2 



E„ 



T 



(24) 



Since E m - m > m v , this is always less than exp[— 0.35(m u /T l/ ) 2 ] for interesting masses m u > T v \ 
unfortunately, the time now appears in the expression for i? m i n , so the dt integral is no longer 
trivial. First, then, let us consider the suppression factor if we ignore the effect of decays 
inside the supernova envelope, integrating from E min = max(m v , k). Then the time integral 
can be done as in the low-mass case, and we can simply write down the time-independent 
suppression factor s = j(m u /T)g(E min /T) ~ j(m v /T)g(m u /T) (if we first set h(k/T) = 1 
when numerically calculating the limits as above). 

These mass effects enable us to break the degeneracy between m v and r, at the price of 
requiring the more complicated analysis of a three-dimensional parameter space. To simplify 
matters, we will base the results for massive neutrinos directly on the limits from the low- 
mass case. That is, we will calculate the limits as before, and then apply the suppression 
factor at the end. We can do this because the suppression factor comes into the expression 
for the flux in exactly the same way as the branching ratio B 1 , so we translate limits on 5 7 
in the massless case to limits on _B 7 x s, where s is the /c-independent part of the suppression 
factor. In addition, we do the calculation for an isotropic decay, and assuming k < m u . For 
two-body decay, this results in the limit 



s5 7 < 3 x 10~ 7 



sB, < 6 x 10- 14 -^^ 
keV sec 



m u r < 10 6 keV sec 
m u T > 10 6 keV sec. 



(25) 
(26) 
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If we allow the effect of decays inside the progenitor envelope, the calculation is somewhat 
more complicated, and the integral over each time bin can no longer be done in closed form. 
We must now recompute everything at each pair of m v and r. We show the results of such a 
calculation for several values of the neutrino mass in Figure || the limits are not too different 
from those with the simpler time-independent suppression. 

For the bremsstrahlung process, the suppression of a high-mass neutrino flux is simpler 
to calculate because the required integral is simply / dE EL#(E), the total energy in the 
massive neutrino species (cf. Eq. pl| we have again assumed k <C T, so the initial integral 
over fM simplifies). Again, we integrate from the same E min = max (A;, E env ); the suppression 
is given by the Boltzmann factor j(E m[n /T). For much of parameter space, this is simply 
the expected j(m u /T). Again, the time-dependence of E env does not change the limits 
significantly. The allowed parameter space for the bremsstrahlung process for a mass of 30 
MeV is 

Tfl T J3 

. v * > 1.5 x 10 12 keV sec or j(mJT)^- < 7 x 10~ 13 key 1 sec" 1 . (27) 
j{m u /T) m v T 



VI. NEARLY DEGENERATE NEUTRINOS 

Thus far, we have assumed that the daughter neutrino in the v —>■ u'j channel is much 
less massive than the parent neutrino. If, however, the mass of the daughter is appreciable, 
the energy of the photon will be decreased by a factor 8m 2 jm 2 = (m\ — m^)/m\. This may 
improve our limits: It would shift the bulk of the photons down from energies too high to 
detect into one or more of the energy channels of the PVO detector. For massless daughter 
neutrinos, of order 1/10 of the photons can be detected; therefore we might expect limits as 
much as an order of magnitude stronger. In fact, for this case, the lower energy window of 
the PVO detector, down to 0.1 MeV, compared with the SMM window, sensitive only above 
4.1 MeV, is actually an advantage. 

To make the matter more precise, we see that in the case of nearly degenerate neutrinos, 
we make the change 

f{E, k)dk -> f[E, {m 2 /5m 2 )k} {m 2 /5m 2 )dk (28) 

where we now are contrained to have photon energies k < (5m 2 /m)E. This, in turn, results 
in changing h{k/T) — > (m 2 /Sm 2 )h[(m 2 /Sm 2 )k/T] « (m 2 /Sm 2 )h(k/T), up to a constant of 
order one. As expected, the flux is enhanced by a factor of (m 2 /8m 2 ) at most energies, so 
we can make the substitution S 7 — > Bj(m 2 /Sm 2 ) in our limits. 

In particular, we are interested in two regimes of Sm 2 /m 2 . To account for the recent 
evidence of a neutrino mass eigenstate, we set m = 3 eV [T(J. The most attractive solution 



to the solar neutrino problem requires Sm 2 = 3 x 10 5 eV 2 , so Sm 2 /m 2 = 3 x 10 6 [11 
To account for the atmospheric deficit, 5m 2 = 1 x 10 _2 eV 2 , so 5m 2 jm 2 = 1 x 10 



|l~2]| . Roughly, we see that our limits will become stronger by factors order 3 x 10 5 and 10 3 , 



respectively. In Figure ^j, we show the limits on the neutrino parameters for these two cases. 
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VII. LONG LIFETIMES 



For long lifetimes (such that the average decay time of the neutrino is comparable to or 
longer than the travel time to the detector), the above formalism becomes too cumbersome, 
because we must integrate over a complicated set of possible paths for the neutrino and 
daughter photon. In this case, we will make several simplifications. At first, we will only 
concern ourselves with the total gamma-ray fluence from the decays, integrated over time. 
Then, we will integrate over the decay time, < < D: 



dN 5 7 



dk 4ttD 2 



J dE L#(E)f(E, k) [l - e - Dm » /ET ] (29) 



where, as before, f(E, k) gives the fraction of neutrinos with energy E decaying into photons 
with energy k. This expression is to be compared with those presented in 0. The cost of the 
simplicity of this expression is the inability to determine the exact time of a photon's arrival. 
For neutrinos and photons travelling on a straight path (/i = 1, appropriate for relativistic 
particles), the arrival time after the supernova light-pulse is £ = £^(1 — v) ~ {t d /2)m 2 /E 2 . 
For long lifetimes, we will be concerned with neutrinos that decay late in their flight: t d ~ D. 
Using this as a typical e-folding time, we have the ansatz that (dN/dt) oc exp (—2tE 2 /Dm 2 ) . 
For £ > (D/2)m 2 /E 2 , this should express the character of the time-dependence. Two effects 
are explicitly missing from this formula: the extra time delay from non-straight paths (of 
the same order as the delay already considered) and the photon energy dependence of the 
time-delay. Moreover, the time-dependence will not have exactly this shape; for short times 
it does not contain the expected slow rise from zero flux, so it is probably safest to use this 
formula integrated over the entire duration of the experiment, and not rely on the detailed 
time evolution, leaving us finally with 



dN 
~dk 



J dE L#{E)f{E, k) (l - e -™&/D"P) (l - e -Dm v /Er^ _ ^ 



Unfortunately, the integration over neutrino energy E is considerably more complicated than 
before, but we can approximate the two exponential decays for various regimes: 

1 _ e -28tEVDml^ ^ _ e -Dm v /Er^ (31) 



1; m < EyJ26t/D, m u /r > E/D 

25tE 2 / Dm 2 ; m > E^26t/D, m u /r > E/D 
Dm u /Er; m < E^26t/D, m v /r < E/D 
2StE/m u r; m > eJqM/D, m v /r < E/D. 



Numerically, these breaks occur at 



I. ( Si \ 1/2 . . 6 . . 1. 



m ~ EJ25t/D ~ 680 eV — - — — ; m v /r ~ E/D ~2x l(T b eV/sec : 



12 MeV \ 8500 sec / ' ' 12 Mev 

(32) 
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where E = 12 MeV is a typical energy for a T v = 8 MeV blackbody. In terms of the 
lifetime r, the latter limit occurs at r ~ Dm u /E « 5 x 10 5 sec (m„/eV) — for masses ~ 1 
MeV, this is roughly r ~ D. When this is proportional to <5t, the flux is approximately 
constant; otherwise, the entire pulse is detected (and its shape is irrelevent). Putting all of 
this together, and doing the integral over E gives 

^h (k/T) 



dN 
dk 



st 



£ 7 120 
AtiD 2 7tt 4 



x < 



mi D 



EtTTLv D 



(33) 



h- X {k/T) 



T 

2^--h,WT) 
T v m u t 

where h n {y) = x n+1 /(l + expx) is similar to h(y) above, and the four cases correspond 
to those in Eq. (|3l]). Here, we have assumed an isotropic distribution of decays in the 
rest frame. As before, these expressions hold for m u < T u and must be modified with the 
appropriate suppression factor otherwise. 

Now, we can just put these results through our statistical machinery and find limits on 
the parameters. We will write the flux as 



dN 
dk 



120 E T 



47rL> 2 7tt 4 T? 



h(k/T) x B~A(m> 



(34) 



where A is the appropriate dimensionless combination of m u and r, along with D, T v and 
5t; the data give us limits on A in each (m u , r) regime (assuming that we can write our 
prior information as a simple probability distribution for A). This gives an approximate 
99% confidence limit of B^A < 1 x 10~ 6 for T v = 8 MeV or, 



B <lx 10~ 6 



8 MeV/ 



m. 



> 155 keVB} /2 



< 1.4 x lO" 15 keV sec" 1 

r 



T 

J- 7/ 



m 



u t>1Ax 10 14 keVsec^ 



8 MeV 

_ZL.V 

8 MeV/ 



m v < 0.4 keV, m v /r > 1.2 x 10~ 9 keV/ sec; 
m v > 0.4 keV,m u /r < 1.2 x 10~ 9 keV/ sec; 

m v < 0.4 keV, m u /r > 1.2 x 10" 9 keV/ sec; 

m v > 0.4 keV, m u /r < 1.2 x 10" 9 keV/ sec. (35) 



Where the regimes overlap, these limits are comparable to those calculated with the more 
detailed models above — because we can only calculate limits on parameters, the details 
of the data and the analysis are unimportant (in fact, the limits of Eq. (^) are stronger 
than, for example, Eq. flT3"|) above; the earlier, more detailed calculation is probably the 
more appropriate limit). Again, for neutrinos with m v T v , these limits are modified with 



B, 



sB^. 



For the bremsstrahlung channel, the flux is changed due to the different kinematics of 
the decay {i.e., the rest frame spectrum of Eq. (|T6|)): 
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dN 
dk 



_ a 2Tl dN 



2-body 



(36) 



(in addition, the functions h n should also be modified to h n+2 ). This is a significant increase 
in flux at for km u < T 2 . As before, we see that the bremsstrahlung spectrum at the detector 
is proportional to 1/k. Now, the limits correspond to B^AT/m < 3 x 10~ 5 or 



m. 



> 2.7 x 10 8 keVR 



T 



m, 



> 800 keVR 1 / 3 



8 MeV 
T 

8 MeV 



m u < 0.4 keV, m u /r > 1.2 x 10~ 9 keV/ sec; 
> 0.4 keV, m u /r < 1.2 x 10" 9 keV/ sec; 



r/S 7 = r e > 1.9 x 10 16 sec m v < 0.4 keV, m u /r > 1.2 x 10~ 9 keV/ sec; 



m, 



> 1.4 x 10 8 keV 



't/B, 



sec 



-1/2 



T 



8 MeV 



^ m v > 0.4 keV, m^/r < 1.2 x 10" 9 keV/ sec. (37) 



VIII. DISCUSSION 

SN1987A not only confirmed astrophysicists' standard model of the Type II (core col- 
lapse) supernovae, but also provided a laboratory for studying the properties of neutrinos. 
The fluence of neutrinos at earth was approximately 10 11 cm -2 per species. This large flu- 
ence and the space-borne gamma-ray detectors operating on SMM and PVO have allowed 
stringent limits to be placed on the radiative decay of neutrinos. 

Although there is only 232 seconds of data, in a single time bin, the SMM detectors 
were sensitive up to energies of 25 MeV, and so would have more likely to detect gamma 
rays from decaying SN1987A neutrinos with a temperature of 4-8 MeV. Because the limits 
are determined by the region of parameter where the background becomes comparable to 
the signal (see the discussion in Sec. |T| above), the long time base and greater resolution 
is actually of little use in determining the limits on the parameters. To show this, we have 
performed our analysis with the SMM data as presented in H , using crude estimates of the 
response matrix and effective area of the detector. As expected, the results are comparable 
to, although slightly weaker than, those quoted in that reference. However, in the case 
where the gamma rays are produced by bremsstrahlung or the neutrino mass states are 
nearly degenerate, the PVO limits are comparable or more stringent. Finally, the amount 
of and quality of the PVO data adds additional confidence to the SMM-based limits. 
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FIGURES 
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FIG. 1. The PVO GBD data for the time immediately before and after the arrival of the light 
from SN1987A at the PVO spacecraft (UT = 27325). Time bins are either 12 or 16 sec; we show 
the average counts per second in each bin, for each energy channel, as marked. 





10 


| 1 1 1 1 | 1 1 1 1 | 1 1 1 1 | 1 

- B 7 =l 


>>,,>>,>!> 




8 














so 








o 


6 


- B 7 =10- 3 \ 






A 


- B 7 =10" 10 \ 
i , , , , i , , , , i , , , A 


,\ , 1 , , i , 1 , 



-5 5 10 15 

log[m 1> T/(keV sec)] 



FIG. 2. The value of the x 2 statistic, defined in the text, as a function of the parameter m u T, 
for values of the branching ratio B 1 as marked, for the decay process v — > 1/7. 
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FIG. 3. Allowed region of the m u T-B^ plane, for the 2-body decay process v — > 1/7, corre- 
sponding to Ax 2 < 2357 (see text). Here and below, contours continue to infinity as long as the 
appropriate assumptions, discussed in the text, still hold. 
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log[m^T e / kev sec] 

FIG. 4. The value of the \ 2 statistic, defined in the text, as a function of the parameter 
777-1, r e = m u r/B e , for the bremsstrahlung process. Also shown is the location of the Ax 2 < 2357 
limit. 
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FIG. 6. Allowed region of the m„T-B 7 plane, for neutrinos of mass 20 MeV, 30 MeV, and 
40 MeV (from left to right), for the 2-body decay process, corresponding to Ax 2 < 2357 (see text). 
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